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Abstract 

Background: Nucleotide pyrophosphatase/phosphodiesterase 7 (NPP7) is the only member of the mammalian 
NPP enzyme family that has been confirmed to act as a sphingomyelinase, hydrolyzing sphingomyelin (SM) to 
form phosphocholine and ceramide. NPP7 additionally hydrolyzes lysophosphatidylcholine (LPC), a substrate 
preference shared with the NPP2/autotaxin(ATX) and NPP6 mammalian family members. This study utilizes a 
synergistic combination of molecular modeling validated by experimental site-directed mutagenesis to explore the 
molecular basis for the unique ability of NPP7 to hydrolyze SM. 

Results: The catalytic function of NPP7 against SM, LPC, platelet activating factor (PAF) and para- 
nitrophenylphosphorylcholine (pNPPC) is impaired in the F275A mutant relative to wild type NPP7, but different 
impacts are noted for mutations at other sites. These results are consistent with a previously described role of F275 
to interact with the choline headgroup, where all substrates share a common functionality. The L107F mutation 
showed enhanced hydrolysis of LPC, PAF and pNPPC but reduced hydrolysis of SM. Modeling suggests this 
difference can be explained by the gain of cation-pi interactions with the choline headgroups of all four substrates, 
opposed by increased steric crowding against the sphingoid tail of SM. Modeling also revealed that the long and 
flexible hydrophobic tails of substrates exhibit considerable dynamic flexibility in the binding pocket, reducing the 
entropic penalty that might otherwise be incurred upon substrate binding. 

Conclusions: Substrate recognition by NPP7 includes several important contributions, ranging from cation-pi 
interactions between F275 and the choline headgroup of all substrates, to tail-group binding pockets that 
accommodate the inherent flexibility of the lipid hydrophobic tails. Two contributions to the unique ability of NPP7 to 
hydrolyze SM were identified. First, the second hydrophobic tail of SM occupies a second hydrophobic binding pocket. 
Second, the leucine residue present at position 107 contrasts with a conserved phenylalanine in NPP enzymes that do 
not utilize SM as a substrate, consistent with the observed reduction in SM hydrolysis by the NPP7-L107F mutant. 



Background 

The nucleotide pyrophosphatase/phosphodiesterase 
(NPP) enzymes are classified as pyrophosphatases or 
phosphodiesterases depending on the type of substrates 
they hydrolyze. NPP1-3 are pyrophosphatases that 
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cleave inorganic phosphates from nucleotides and their 
derivatives [1]. In contrast, NPP2, NPP6 and NPP7 are 
phosphodiesterases that hydrolyze phosphodiester bonds 
in lipids and their derivatives [2,3]. NPP4 and NPP5 are 
yet to be characterized in terms of substrates and types 
of activity. NPP2 is the only one of the seven mamma- 
lian NPP family members that is known to possess both 
pyrophosphatase and phosphodiesterase activities [4]. 
NPP2 is also the only mammalian family member to 
have been characterized by X-ray crystallography [5,6]. 
As a phosphodiesterase it exhibits a lysophospholipase 
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D (LPLD) activity, in contrast to NPP6 and NPP7 which 
exhibit lysophospholipase C (LPLC) activity. Besides 
these specific activities, NPP family members also exhi- 
bit distinct substrate preference profiles but no clear 
structural basis for the preferences has been proposed 
[1]. NPP family members are of considerable interest 
due to their role in biological functions ranging from 
bone mineralization to cancer [4,7,8]. NPP7 specifically 
has been linked with anti-inflammatory and anti-tumori- 
genic activities through its influence on the conversion 
of SM to ceramide [9-11], however homozygous knock- 
out of NPP7 in mice did not lead to spontaneous 
tumorigenesis although lipid digestion and absorption 
was altered [12]. Mutagenic deletion of a segment from 
NPP7 including H353, a residue involved in chelating 
one of the essential divalent metal cations, has been 
identified in the rapidly proliferating colon cancer HT- 
29 cell line [13]. 

NPP7 hydrolyzes lysophosphatidylcholine (LPC, Figure 
1), platelet activating factor (PAF, Figure 1) and sphin- 
gomyelin (SM, Figure 1) with a lysophospholipase C 
activity but exhibits a preference for SM over LPC and 
PAF [10]. Several studies have examined the roles of 
specific amino acids in the hydrolysis of SM. Three of 
the conserved metal-chelating residues from the NPP 
family (D199 and D246 [10,14], as well as H353 [13]) 
have been mutated to alanine in NPP7, and each mutant 
was devoid of sphingomyelinase activity. The conserved 
catalytic threonine residue, T75, has also been mutated 
to alanine with the expected loss of sphingomyelinase 
activity [14]. Additional amino acids near the catalytic 



residue have been mutated in attempts to change sub- 
strate specificity to include nucleotide pyrophosphates, 
including M74K, S76F, and C78N [14]. These mutations 
all eliminated sphingomyelinase activity without produ- 
cing pyrophosphatase activity. It is unclear whether the 
mutated proteins were properly folded, although expres- 
sion levels were comparable to wild type. More recently, 
mutations based on a comparative model of NPP7 were 
performed [15]. In contrast to the result obtained upon 
mutation of M74 to K, mutation of M74 to L enhanced 
sphingomyelinase activity, again without gain of pyro- 
phosphatase activity. The model displayed close interac- 
tion between F275 and the choline methyl groups, best 
described as a cation-n interaction. Mutation of F275 to 
glycine nearly eliminated sphingomyelinase activity. The 
model displayed no close interaction between F141 and 
substrates, however, the F141S mutation showed activity 
nearly as poor as the F275G mutant, indicating either 
incorrect selection of the conformation of the loop 
involving residues 140-171, or incorrect substrate posi- 
tioning. Therefore the interactions between NPP7 and 
SM are partially defined, and modeled interactions with 
other substrates have not been supported by experimen- 
tal results. Further efforts are clearly needed to develop 
a complete and coherent picture of substrate recognition 
and discrimination by NPP7. 

Lyso-platelet activating factor (lyso PAF), sphingosyl- 
phosphorylcholine (SPC), and para-nitrophenyphosphor- 
ylcholine (pNPPC, Figure 1) are high affinity substrates 
for NPP7 (Tables 1, 2, 3 and 4, and [10]). Here we 
examine the amino acid residues within the NPP7 active 




LPC 16:0 



PAF 16:0 



SM 16:0 



pNPPC 



Figure 1 NPP7 substrate structures. 
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Table 1 Kinetic parameters for hydrolysis of LPC 16:0. 



Enzyme 


K m (MM) 


V max (MM/sec) 
(x 10" 2 ) 


kcattsec 1 ) 


UKm (MW 1 ) (x 10 4 ) 


1 1 07F 
L I U/ r 


^1 J. TO 
O 1 X z.o 


j.J X U.U3 


^-.Z JZ U.U3 


a. o _i_ i a. 
O.o jz I .D 


Y142A 


69 ± 0.9 


3.0 ± 0.01 


3.6 ± 0.01 


5.2 ± 0.05 


Wild type 


57 ± 0.7 


2.5 ± 0.01 


2.9 ± 0.01 


5.1 ± 0.1 


E169Q 


88 ± 6.8 


3.3 ± 0.02 


3.8 ± 0.03 


4.3 ± 0.4 


E169A 


100 ± 15 


1 .0 ± 0.04 


1 .2 ± 0.05 


1 .2 ± 0.2 


F80A 


120 ± 3.0 


0.7 ± 0.01 


1.4 ± 0.01 


1 .2 ± 0.3 


Y166A 


100 ± 5.0 


0.9 ± 0.03 


1.1 ± 0.03 


1 .0 ± 0.2 


F275A 


No activity 









site that may be involved in substrate recognition by 
NPP7. We employed homology modeling to predict pro- 
tein-substrate interaction points followed by experimen- 
tal validation of the predictions by mutagenesis and 
biochemical assays. We used the bacterial Xac NPP 
crystal structure [16] as well as two recently reported 
NPP2 crystal structures [5,6] to generate homology 
models of NPP7 for this work. Our modeling results 
reveal a common choline headgroup binding pocket for 
all substrates. The strong interactions between the catio- 
nic choline and the tt electrons of the aromatic side- 
chain of F275 were easily identified by standard docking 
algorithms. Experimental mutagenesis results are consis- 
tent with the predicted choline headgroup positions for 
all substrates. The hydrophobic tails proved more chal- 
lenging to place in the complex. Docking algorithms 
typically cannot exhaustively search the conformations 
and positions of long and flexible molecules lacking any 
polar groups that would be involved in directional inter- 
actions such as hydrogen bonding. Manual placement 
was used, and refinements of hydrophobic tail positions 
by molecular dynamics were required in order to corre- 
late the substrate complex models with the experimental 
mutagenesis results in the hydrophobic tail region. The 
molecular dynamics simulations were important to 
demonstrate the high mobility of the hydrophobic tail 
positions, as well as to refine the initial complex of 



sphingomyelin (SM) to explain the different pattern of 
experimental mutation impacts on activity for this sub- 
strate, relative to the other three examined. 

Methods 

Molecular modeling 

Molecular modeling and simulation was done using the 
Molecular Operating Environment (MOE) software 
package [17]. Three homology models of NPP7 were 
built in MOE by aligning the human protein sequence 
[2] (Genbank code AY20633) with either the bacterial 
Xac crystal structure from the protein data bank (PDB 
[18] entry 2GSO[16]), the rat NPP2 crystal structure 
(PDB entry 2XR9[5]) or the mouse NPP2 crystal struc- 
ture (PDB entry 3NKM[6]). The homology models were 
energy minimized using the MMFF94x force field, a 
derivative of the MMFF94 [19] force field with improved 
treatment of planar nitrogens. Minimizations were ter- 
minated when the energy gradient reached a root mean 
square (RMSG) value of 1 kcal/mol-A in order to pre- 
vent collapse of binding pockets to enhance hydropho- 
bic interactions in the absence of ligand. 

Substrate docking with the MOE software was used to 
study interactions between NPP7 and its previously 
known substrates and their closely related analogs. 
Docking used placements of ligand atom triangles into 
protein site alpha sphere triangles with subsequent 



Table 2 Kinetic parameters for hydrolysis of PAF 16:0. 



Enzyme 


K m (MM) 


V max (|4M/Sec) 

(x 10" 2 ) 


k cat (sec 1 ) 


k cat /K m flvrW 1 ) (x 10 4 ) 


Y142A 


75 ± 2 


2.0 ± 0.03 


2.5 ± 0.04 


3.3 ± 0.6 


L107F 


94 ± 3 


2.4 ± 0.3 


2.9 ± 0.3 


3.1 ± 0.2 


Wild type 


82 ± 2 


1.3 ± 0.02 


1.6 ± 0.02 


2.0 ± 0.8 


E169Q 


130 ± 8 


1.9 ± 0.1 


2.3 ± 0.2 


1.7 ± 0.2 


F80A 


120 ± 3 


1.2 ± 0.03 


1.4 ± 0.03 


1.2 ± 0.4 


E169A 


95 ± 6 


0.8 ± 0.02 


0.80 ± 0.02 


1.0 ± 0.2 



Y166A 200 ±10 0.7 ± 0.02 0.80 ± 0.04 0.4 ± 0.2 



F275A 



No activity 
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Table 3 Kinetic parameters for hydrolysis of SM. 



Enzyme 


K m (MM) 


V max (MM/sec) 
(x 10" 2 ) 


kcat (sec 1 ) 


k cat /K m (MW 1 ) (x 10 4 ) 


pi Aon 


1 ^ + ^ 

I J X J 


D.Z JZ U.Uo 


D.O JZ U. I 


ZJ.U 2Z D.U 


Y166A 


12 ± 0.4 


2.6 ± 0.01 


3.1 ± 0.2 


25.8 ± 5.0 


F80A 


14 + 6 


2.2 ± 0.3 


2.6 ± 0.3 


18.6 ± 5.0 


Wild type 


23 ± 3 


3.0 ± 0.1 


3.7 ± 0.2 


16.5 ± 4.0 


E169A 


17 + 2 


1.8 ± 0.1 


2.2 ± 0.2 


11.6 ± 2.0 


L107F 


37 ± 3 


3.5 ± 0.1 


4.2 ± 0.1 


11.4 ± 4.0 


Y142A 


47 ± 7 


3.1 ± 0.2 


3.7 ± 0.2 


7.9 ± 3.0 


F275A 


12 ± 1 


0.7 ± 0.1 


0.9 ± 0.2 


7.5 ± 2.0 



focefield refinement using MMFF94x [19]. One hundred 
initial placements were used in each docking run. 
pNPPC was directly docked into the active site as a 
whole molecule since it is small and has limited flexibil- 
ity. The phospholipid substrates; LPC, PAF, and SM, are 
highly flexible because of long hydrocarbon chains in 
their structures. For this reason stepwise docking was 
done starting with the polar head groups of each sub- 
strate. Wall restraints were used to keep the substrate 
phosphate within 8 A of the enzyme's active residue, 
threonine 75. The hydrophobic chains were then 
extended three carbon atoms at a time followed by 
energy minimization for the newly added atoms at each 
step. The protein surface feature within MOE was used 
to visualize the various channels that could accommo- 
date the substrate hydrophobic tails during binding. 
Bond rotations were done on the extended chains fol- 
lowed by energy minimization to explore the best posi- 
tions for each chain and thus the minimum energy 
conformation for each substrate. The lowest energy con- 
formation was then selected for each substrate. 

Molecular dynamics was performed on selected NPP7 
substrate complexes using the MOE software. Simula- 
tions were performed using the MMFF94x forcefield 
with the Generalized Born implicit solvation model [20]. 
The NVT ensemble was used with constraints on bonds 
to hydrogens. Atomic positions were recalculated at 2 fs 



timesteps using the Nose-Poincaire-Andersen equations 
of motion [21,22]. A heating phase of 100 ps was used 
to raise the simulation temperature from 0 K to 300 K. 
The production phase began at the end of the heating 
phase and continued to the 2000 ps timepoint. Snap- 
shots were saved every 1 ps during both phases of the 
simulation for analysis. 

NPP7 Plasmid design 

A mammalian expression vector pcDNA4/TO/myc His 
B containing the full-length human NPP7 sequence was 
a generous gift from Dr. Rui-Dong Duan (Lund Univer- 
sity, Sweden). A FLAG -tag and a stop codon were incor- 
porated at position 415 using the sense primer 5'-CCC 
ATG CTG CAC ACA GAC TAC AAG GAC GAC 
GAT GAC AAGTAGGAATCTGCTCTTCCG-3' and 
antisense primer 5'-CGGAAG AGCAGATTCC- 
TACTTGTCATCGTCGTCCTTGTAGTCTGTGTG- 
CAGCATGGG-3'. This generated NPP7ex-FLAG which 
had previously been demonstrated [23] to retain cataly- 
tic activity and localize in the cell culture medium. 
Mutagenesis was done by PCR using the Quickchange® 
site directed mutagenesis kit from Qiagen-USA. This 
customized gene was further sub cloned into the mam- 
malian vector, pcDNA3.1 (+) at the NotI and BamHI 
restriction sites. All mutants were generated from this 
customized plasmid using Qiagen's multi site-directed 



Table 4 Kinetic parameters for hydrolysis of pNPPC. 



Enzyme 


K m (MM) 


V max (MM/sec) 
(x 10' 2 ) 


kcat (sec 1 ) 


kcat/K m (NT 1 sec 1 ) (x 10 4 ) 


L107F 


310 ± 5.0 


1 .4 ± 0.05 


1.7 ± 0.1 


0.60 ± 0.02 


Wild type 


360 ± 8.0 


1.1 ± 0.01 


1.4 ± 0.1 


0.40 ± 0.03 


Y142A 


400 ± 4.0 


1.4 ± 0.01 


1.7 ± 0.1 


0.40 ± 0.02 


F80A 


450 ± 10.0 


0.63 ± 0.04 


0.8 ± 0.1 


0.20 ± 0.1 


E169A 


170 ± 30.0 


0.42 ± 0.03 


0.5 ± 0.03 


0.30 ± 0.02 


E169Q 


440 ±5.0 


1 .2 ± 0.02 


1.4 ± 0.1 


0.30 ± 0.02 


Y166A 


820 ± 20.0 


0.62 ± 0.03 


0.7 ± 0.04 


0.09 ± 0.01 


F275A 


120 ± 5.0 


0.23 ± 0.03 


0.3 ± 0.03 


0.03 ± 0.01 
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Table 5 Sense primers for hNPP7ex-FLAG mutations. 



MDD"7 

INrr/ 


Sense primer 


mutation 




F80A 


TGACCAGCC CTGCCACGCC ACCCTGGTCA 




CCGGCAAATA 


L107F 


ACCACCAGCA AGGTGAAG TTC 




CCCTACCACGCCACGCTGG 


Y142A 


TGGCTCCTO TCGCCCCGG GCGGGAACGT 


Y166A 


GGCATCGCAC ACAACGCCAAAAATGAGACG 


E169A 


TCGCAC ACAACTACAAAAATGCGACG GAGTGGAGAG 


E169Q 


TCGCAC ACAACTACAAAAATGAGACG GAGTGGAGAG 


Y194A 


ATCTGGTCACACTCGCCTTCGGGGAGCCGGACTCC 


F275A 


CGGGACATCG AGmGAGCT CCTGGACTAC 



mutagenesis protocol using only the sense primer. All 
sense primers used in generating our various mutations 
are shown in Table 5. All PCR products were amplified 
in Nova blue E. coli cells and purified according to the 
manufacturer's protocols (EMD Chemicals). 

NPP7ex-FLAG expression 

Human embryonic kidney cells (HEK293T) were seeded 
in Dulbecco's modified Eagle's medium (DMEM) con- 
taining 10% heat inactivated fetal calf serum and 2 mM 
L-glutamine. The cells were grown overnight at 37°C 
and 5% C0 2 up to a 80% confluence before transfection 
with purified hNPP7ex-FLAG insert subcloned in the 
pcDNA3.1(+) mammalian vector. Transfection was done 
in the presence of polyfect transfection reagent from 
Qiagen according to the manufacturer's protocol. Five 
to eight hours after trasfection, the culture medium was 
changed to serum free DMEM and cells incubated for 
48 hours at 37°C and 5% C0 2 . Conditioned media con- 
taining secreted hNPP7ex-FLAG protein was collected 
and concentrated using 10,000 molecular weight cutoff 
(MWCO) filters (Millipore, Billerica, MA, USA) by cen- 
trifugation at 3000 x g for 20 minutes. Concentrated 
NPP7 protein was stored at 4°C in NPP7 assay buffer, 
pH 8 without addition of adjuvants. 

Western blots 

Western blots with the M2 Anti-FLAG antibody (Sigma- 
Aldrich, St. Louis, MO, USA) were used to verify 
expression of hNPP7ex-FLAG and its mutants by 
HEK293T cells. Following Tris-HCl-SDS-PAGE (4-15%) 
proteins were transferred onto poly-m-vinylidene 
difluoride (PVDF) membrane (Bio-Rad, Hercules, CA, 
USA) by electroblotting at 100 V in Tris/Glycine SDS 
buffer for 1 hour. Nonspecific binding sites were 
blocked in 50 mM Tris buffered saline (TBST: 2.7 mM 
KC1, 0.138 M NaCl, 0.5% (v/v) Tween®20) and 5% (w/v) 
non-fat dry milk (NFDM). The membrane was then 
treated with Anti-FLAG M2 antibody (Sigma-Aldrich, 



USA) (1:2000 in TBST + 3% NFDM) and anti-goat IgG 
antibody (Sigma) (1:5000 in TBST + 3% NFDM). Finally, 
peroxidase-labeled secondary antibodies were used for 
chemiluminescence detection. 

Normalization and quantification of hNPP7ex-FLAG in 
CCM 

All NPP7 activity assays were performed with concen- 
trated conditioned media (CCM) quantified at least 
twice by western blot analysis. First, a scanning western 
blot was run with varying dilutions of the CCM immedi- 
ately after concentration to assess differences in the 
amounts of protein in the CCM. Based on the results of 
the first western blot, all samples were diluted to the 
same concentration then again verified by a second 
immunoblot. A known concentration of the carboxy- 
terminal FLAG-BAP fusion protein from Sigma- Aldrich- 
USA was included on the second western blot as a stan- 
dard to quantify the amounts of protein expressed for 
each sample. All protein bands were normalized relative 
to the FLAG -BAP standard using the TotalLab software 
program on a Fotodyne imager (Hartland, WI). 

Assays of NPP7 catalytic activity 

All assays of NPP7 catalytic activity were performed 
three independent times with triplicate samples in each 
experiment. Standard deviations of results from three 
independent experiments provide error estimates. 
Fluorescence-based assay 

Hydrolysis of all three lipid substrates by NPP7 was 
monitored using a modification of Invitrogens Amplex® 
Red phosphatidylcholine-specific phospholipase C assay 
kit. NPP7 assay buffer was prepared as previously 
described [2] containing 50 mM tris-HCl, pH 8.0, 150 
mM NaCl, and 10 mM taurocholic acid but without 2 
mM EDTA. Each well on the plate was loaded with a 
total of 60 uL consisting of 20 uL each of substrate, 
Amplex red cocktail (10 uM Amplex Red reagent, 0.1 
U/mL choline oxidase, 1 U/mL horseradish peroxidase, 
4 U/mL alkaline phosphatase) and normalized concen- 
trations of enzyme diluted from CCM. NPP7 activity 
was monitored by the level of fluorescence using excita- 
tion and emission wavelengths of 571 nm and 585 nm, 
respectively. Data were collected at 1 minute intervals 
over a 2 hour period in a 96-well, half area plate (Corn- 
ing Inc., Corning, NY) using a Synergy 2 multi-well 
plate reader (BioTek, Winooski, VT, USA). 
Absorbance-based assay 

Hydrolysis of pNPPC was monitored using an absor- 
bance assay system based on the ^m-nitrophenolate 
product that has a maximum absorbance at 405 nm. 
Each well on the plate was loaded with a total of 60 uL 
consisting of 20 uL each of NPP7 assay buffer, substrate 
and normalized amounts of enzyme diluted from CCM. 
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NPP7 activity was monitored by the level of absorbance 
at 405 nm every minute for a 2 hour period in a 96- 
well, half area plate (Corning Inc., Corning, NY) using a 
Synergy 2 multi-well plate reader (BioTek, Winooski, 
VT, USA). 

Calculation of kinetic parameters 

Relative fluorescence and absorbance units were con- 
verted into product concentrations. For the fluorescence 
and absorbance assays, varying concentrations of resoru- 
fin or para-nitrophenol were used to generate standard 
curves. The trend line equation from a plot of standard 
concentration versus fluorescence/absorbance was used 
to convert all relative fluorescence and absorbance units 
to product concentrations. Curve fitting was done with 
the kaleidograph software program to obtain K m and 
V max - K ca t was calculated using the V max obtained and 
8.33 nM enzyme concentration used in all assays. 

Results 

hNPP7 homology models 

Three homology models of hNPP7 were developed 
using different crystallized NPP structures. Crystal 
structures of a bacterial NPP were first reported in 
2006 [16]. This year, additional crystallographic struc- 
tures of NPP2 from two different rodent species 
became available [5,6]. Notably, the bacterial NPP and 
mouse NPP2 structures have been crystallized in the 
presence and absence of products, AMP and LPA, 
respectively. Figure 2 shows that there is essentially no 
induced fit of the enzyme to these product structures. 
The similarity between structures crystallized in the 



presence and absence of product can be quantitated by 
heavy-atom root mean square deviation (RMSD) values 
for atoms within 4.5 A of the product. The RMSD 
values are only 0.165 A for the bacterial NPP and 
0.438 A for the mouse NPP2. Thus either the free or 
product-bound crystallographic structures should serve 
equally well as a homology modeling template. Amino 
acid sequence identity is similar between hNPP7 and 
the three available NPP structures (24.7-25.3%). How- 
ever, alignments against the available templates show 
substantial differences in the number of amino acids 
that align against gaps in the template structures. The 
longest stretch of hNPP7 amino acids aligned against 
gaps in 2GSO [16], 2XR9 [5], and 3NKM [6] were 6, 
14, and 14, respectively (Figure 3). Models developed 
from each template structure are shown in Figure 4, 
emphasizing regions in the model that aligned against 
gaps in the template structure as well as the position 
of N146, a confirmed glycosylation site [14]. The long- 
est gap-aligned stretches all occur in the hNPP7 
sequence between P143 and E185. Due to the much 
smaller gap-aligned segments in the 2GSO-based 
model of hNPP7, this model was used in all docking 
studies with substrates. Two prior modeling studies of 
hNPP7 also used a crystallographic structure of the 
bacterial NPP as the modeling template [15,24], 
although one of the studies was published prior to the 
availability of the NPP2 crystal structures [15]. The 
insertion of residues in the NPP7 sequence within the 
segment interacting with LPA in NPP2 limits the uti- 
lity of the NPP2 crystal structures as templates even 
though they share a common substrate with NPP7. 





Figure 2 Comparison of NPP family member crystal structures in the absence (magenta) and presence (blue) of ligand A Comparison 
of Xac. NPP crystal structures 2GSO and 2GSU (with AMP). All-atom RMSD for residues within 4.5 A of AMP is 0.017 A. B. Comparison of mouse 
NPP2 crystal structures 3NKM and 3NKP (with LPA). All-atom RMSD for residues within 4.5 A of LPA is 0.44 A. Glu308 is labeled to emphasize the 
sidechain showing the largest difference between the free and bound forms of the enzyme. 
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2GSC T PHALLLI SIDG1RADMLDR — GIT FM1SH E^EEGVR&RWM&FS YPSLT FFNH YTL VTGLRFDHHGIVHNSMRDPTL 

2 XK 9 VRP PLI I FS VDGFRAS YMKKGS KVM PN lEKLfl SC GTHAP YMRPV YP TKT FPNL YTLAT GLYPESHG I VGNSMY DP VF 

3NKM VRP PL 1 1 F S VDGFRA S YM KKGS KVM PN I EK LR SC GTHAP YMRPV Y P T KT F P N L YTL AT GL Y PE S HG IV GNS MY D P VF 

h N F F 7 SQNKLLLVSFDG FRW NYD QD — VDT PN LDAMA RDGVKAR YMTPAFV TMT SPCHFTLVTGK YIE NHGWH -NMY YN TT 



*::.*.**:*. **:. : ^ - * # ; * * _ * * ^ * 

2GS0 - GG FW L SKS - EAVGD AR WW — G GE PVH VG VENT G QHAAT ft S ?f P G SE AAI KG VRP SQ W RH YQK GVRLDTR VD A 

2XR9 DAS FHLR- G-REKFNHRWWG— GQPLW IT ATKQGVRAGT FF WSVSIPHERRILT 

3NKM DAT FH LR- G- REK FN HRWWG — GQPLWIT ATKQG VR AGT FF — -WSV5I PHE RRILT 



hNPP7 - SKVKL P Y H - ATLGI QR WWDNG S V P I W I T AQRQG LRAG SFFYPG GN VT YQG VAV TRSRKEGIAHW YKK ET EWRAW IDT 

■# * > - # - ■ •■ * 

2GS0 VRGVf LAT DG AQRMR LVTL YFE HVDE AGH DHG PE SRQ YA DA VRA VDAA I GR LL AGM QR DG TRA RT N 1 1 W SDHGMAE VA 

2XR9 ILQWLSLPDNERPS VY AF YSE QP DF S GHK YG PFGPEMT P LRE X DKT VGQ LMDGL KQLRLHR CVNV I FV GDHGMED VT 

3NKM ILQWLSLPDNERPS VYAF YSEQFDF SGH KYGPFGPEMTNPLRE IDKT VGQLMDGLKQLKLHRCVNV I FV GDHGMED VT 

hNPF7 VMAWFTEE DLDLVTLYFGEPDS TGHRYGPESPE RREMVRQVDRT VGYLRES I ARNHLTDRLNL I IT SDHGMTT VD 

2GS0 PGHAI SV EDI APPQI ATAIT DGQVI GFEP LPG — QQAAAEA 5 VLGAH DH YDCWRKA EL PAR HQ Y GS HPR 

2XR9 CDRTEFLSNYLTNVDDITLVPGTLGRIRA KSIYDPKTI IAALTCKKPD— QKFKPYMKQHL PKRLHYANNRR 

3NKM CDRTE FL S NY LTNVDDI TLVPGTLGRIRPK — IP NNL KYDPKAI I ANL TCKKPD--QHFKPYMKQHL PKRLHYANNRR 

r.NPP7 KRA GD L V E FH K FPN FT FR DI E FE LL D Y G P NG- ML LPKEG RL EKVY DAL KDAH PKLHVYKKEAF PEAFH YANNPR 



2G30 I PSL VCQMH EGWDA L F PDK LAKRA QRG TRG SHG YD PAL P SMRAV FL AQG PDL AQGKT L PG FDNVDV YALMS RL LGI 

2XRS IEDI HLL VDRRWH V AR KP C F FQG DH G FDNK VN SMQT V FVG YG FT FKY RTK VPPFENIE L YN VMC DL LGL 

3NKM I E DL HLL VE RRWHVARKPL DV YK KPS GKC F FQG DH G FDNK VN SMQT V FVG YG FT FKY RT K V P F FEN I E I YN VMC DLLGI 

h>JFF7 VT — PLLMYSDLGYV IHGRIN FNNG EH G FDNK DM DMKT I FRAVG PS FRAGLE VE P FE S VH V YE LMC RL LGI 

; ; m * + *:::* . ** : : 

2GSO PAAP NDG NP AT LLP ALRM 

2XR 9 KP A P NNGTH GS LNH LLRT 

3NKM KPAPNNGTHGS LNH LLRT 

hNPP7 VPEANDGHLATLLPMLHT 
. . * : * - : + * : 

Figure 3 Alignment of NPP enzyme crystal structures and hNPP7. Positions with conserved residue identity are shown in blue and marked 
with asterisks. Positions with conservative substitutions are shown in green and marked with colons. 



hNPP7 substrate complexes 

Three natural NPP7 substrates (PAF 16:0, LPC 16:0, and 
SM 16:0) and one synthetic substrate (pNPPC) were 
docked into the model of hNPP7. The natural substrates 
were flexibly docked with shortened hydrophobic tails 
and then extended manually. The synthetic substrate 
was docked flexibly. The docked positions of each sub- 
strate are shown in Figure 5 and 5a three dimensional 
model of all four in complex with the hNPP7 homology 
model are provided in additional file 1. The four com- 
plexes show a common position of the choline head- 
group near F275 (cation-n interaction), and a common 
phosphate group position interacting strongly with one 
zinc ion (chelated by D199, H203, and H353, shown to 
the left in all figures) and moderately with the second 
(chelated by D39, D246 and H247, shown to the right in 
all figures) (ionic interactions). All four substrates 
extend from the phosphate group along the bottom 
edge of the parallel p -sheet bundle (van der Waals inter- 
actions). Only the short acetyl sidechain of PAF 16:0 
(Figure 1A) and the long palmitoyl chain of SM 16:0 
(Figure ID) extend above the edge of the parallel P- 



sheet bundle (van der Waals interactions). These com- 
plexes were used to select a set of mutation sites to pro- 
vide experimental feedback and validation of the 
models. The closest heavy atom distances from each of 
these selected mutation sites to the docked substrates 
are shown in Table 6. Based on these complexes, muta- 
tion of F275 to a non-aromatic amino acid was expected 
to result in a less active mutant against all four sub- 
strates. L107 was selected for mutation to an aromatic 
residue with the expectation that this could introduce a 
favorable cation-n interaction with the choline group, 
resulting in a more active mutant against all four sub- 
strates. However, L107 also makes van der Waals con- 
tact with the sphingoid chain of SM 16:0 (3.9 A) and 
the more rigid aromatic ring geometry might occlude 
the space occupied by this chain. Three residues located 
at various positions near the common position occupied 
along the bottom edge of the parallel p -sheet bundle 
were selected for mutation, including F80, Y142, and 
Y166. One mutation site was selected near the chains 
extending above the parallel P-sheet bundle, E169. 
These latter four mutations were expected to show 
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Figure 4 Comparison of hNPP7 models based on three different template structures. Green ribbons represent hNPP7 residues aligned 
against residues in the template. Magenta ribbons represent hNPP7 residues aligned against gaps in the template. Zn 2+ ions are shown as blue 
spheres and the position of N 146 is shown as a spacefilling model. A. hNPP7 based on PDB entry 2GSO[16] {Xac NPP), 24.7% amino acid 
identity. B. hNPP7 based on PDB entry 2XR9[5] (Rat NPP2), 24.9% amino acid identity. C. hNPP7 based on PDB entry 3NKM[6] (Mouse NPP2), 
25.3% amino acid identity. 
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Figure 5 Comparison of docked substrate positions. Zn 2+ ions are shown as light blue spheres. Substrates are shown as ball & stick models. 
hNPP7 residues selected for mutagenesis studies are shown as green sticks and labeled. A. PAF 16:0. B. LPC 16:0. C. pNPPC. D. SM. 



variable differences in activity relative to wild type NPP7 
when tested using different substrates due to the vari- 
able contact distances displayed in these preliminary 
docked complexes. For example, Y166 interacted closely 
with the sphingoid tail of SM 16:0 (3.6 A) but showed 
negligible contact with pNPPC (7.4 A). 

hNPP7 mutations 

hNPP7ExFLAG (NPP7-WT) was obtained from the 
medium of transfected HEK-293T cells without further 
purification. To ensure that these cells did not produce 
and export proteins or factors that would interfere with 



Table 6 Closest distances (A) between docked substrates 
and residues selected for mutation in initial docked 
complexes. 





PAF 16:0 
tail 1-16:0 
tail 2-2:0 


pNPPC 


LPC 16:0 


SM 16:0 
tail 1 -sphingoid 
tail 2-16:0 


Choline interactions 


F275 


4.5 


3.7 


3.6 


3.8 


L107 


7.2 


7.5 


5.6 


8.2 


Tail 1 Interactions 


L107 


4.5 


6.2 


5.1 


3.9 


F80 


6.2 


9.5 


6.1 


7.1 


Y166 


4.8 


7.4 


4.8 


3.6 


Y142 


6.9 


6.5 


6.9 


8.8 


Tail 2 Interactions 


E169 


4.0 






5.8 



the Amplex Red assay, fluorescence readings using assay 
buffer (AB), CCM from mock-transfected cells, and 
CCM from NPP7-transfected cells were compared (Fig- 
ure 6). Figure 6 demonstrates that only the CCM from 
NPP7-transfected cells produced a fluorescent signal 
using the Amplex Red assay, and only when alkaline 
phosphatase was present to convert phosphocholine 
produced by LPLC activity to choline. There was no 
fluorescent signal in the absence of alkaline phosphatase, 
eliminating any contribution to the fluorescent signal 
due to direct production of choline by LPLD activity. 
NPP7-WT and selected mutants were then character- 
ized for their hydrolytic activity against PAF 16:0, LP A 
16:0, pNPPC, and SM 16:0. The kinetic parameters 
obtained from these assays are provided in Tables 1, 2, 
3, and 4, and a graphical comparison of relative activity 
is shown in Figure 7. 

Two mutants, F275A and L107F, were selected on the 
basis of their proximity to the choline headgroup of the 
substrates (Table 6). Relative to NPP7-WT activities, the 
F275A mutant showed either no catalytic function (LPA 
16:0 and PAF 16:0), a 13-fold reduction in k cat /K m value 
(pNPPC) or only a 2.2-fold reduction in k cat /K m value 
(SM 16:0). These results are consistent with the models 
which suggested that the F275A mutation should be 
detrimental to the recognition and hydrolysis of all sub- 
strates. In contrast to the consistently reduced activity 
of the F275A mutant, the L107F mutant activity differed 
from WT in a substrate-dependent fashion. The L107F 
mutant hydrolyzed PAF 16:0, pNPPC, and LPA 16:0 
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Figure 6 Fluorescence signals from Amplex Red assays of LPC hydrolysis in the presence (grey bars) and absence (white bars) of 
alkaline phosphatase using assay buffer (AB) or conditioned media from either mock or NPP7-transfected HEK293T cells. 



between 1.3 and 1.6 times more effectively than NPP7- 
WT. Conversely, SM 16:0 was hydrolyzed 1.4 times less 
effectively by L107F relative to NPP7-WT. These experi- 
mental results are consistent with the model-based sug- 
gestion that replacement of the nonpolar isobutyl 
sidechain of leucine with the aromatic sidechain of 



phenylalanine allowed formation of an additional cation- 
tt interaction with the substrates, but produced crowd- 
ing against the SM sphingoid tail. 

Three mutants, F80A, Y166A, and Y142A, were 
selected on the basis of their proximity to a hydrophobic 
tail in all substrates (Table 6). These mutants showed 




L107F Y142A 



E169Q 



E169A 



F80A Y166A F275A 



Mutation 

Figure 7 Relative k cat /K m values of all mutants for hydrolysis of all substrates. Substrates include LPC (white bars), PAF16:0 (light grey bars), 
pNPPC (dark grey bars) and SM (black bars). The relative k cat /K m values were obtained by dividing the k cat /K m value of each mutant with that of 
wild type enzyme against the same substrate. 
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Figure 8 LPC 16:0 position in NPP7 relative to mutated residues. Panel A. Structure shown represents an energy-minimized snapshot from 
the MD simulation. Distances from sidechains to LPC 16:0 are 2.7 A (E169 to hydroxyl), 2.9 A (L107 to choline), 3.0 A (F275 aromatic centroid to 
choline), 3.0 A (F80 to nonpolar tail), 7.2 A (Y166 to nonpolar tail), and 9.0 A (Y142 to nonpolar tail). The two distances showing the greatest 
standard deviations when measured at 1 ps intervals during a 2000 ps simulation trajectory were to F80 (2.0 A) and to Y166 (1.7 A). Panel B. 
Molecular dynamics simulations of the LPC 16:0 complex with NPP7 show dynamic motion of the nonpolar tail within the hydrophobic channel 
as reflected in the distances between the nonpolar tail and the sidechains of F80 and Y166. This dynamic motion results in close interactions at 
any given time point with either F80 or Y166. 



opposing patterns of relative impact on SM relative to 
the remaining three substrates. The Y142A mutant 
showed only half the activity of the WT enzyme against 
SM, but showed similar or better activity against PAF, 
LPC and pNPPC (Tables 1, 2, 3, and 4). The F80A and 
Y166A mutants showed reduced activity against LPC, 
PAF and pNPPC, but showed similar or better activity 
against SM (Tables 1, 2, 3, and 4). This pattern of activ- 
ity differences was not predicted by the initial model 
complexes of these substrates. However, the hydropho- 
bic tails of the endogenous substrates had been placed 
manually, and exploration of alternate tail positions had 
been limited. Molecular dynamics simulations of LPC 
16:0 and SM were therefore used to allow the substrate 
tails to explore nearby low-energy placements within the 
enzyme for a representative substrate from each pattern 
of activity differences. Molecular dynamics simulations 
showed a time-dependent alternating interaction of the 
LPC 16:0 hydrophobic tail with F80 and Y166 based on 
distances that oscillated between 2.5-3.0 A and greater 
than 6.5 A (Figure 8). This result is consistent with sub- 
stantially reduced activity due to mutation at either site 
and suggests that the hydrophobic tail is not restricted 
to a confined location, which would be an entropically 
unfavorable binding contribution. Instead, the hydro- 
phobic tail retains considerable flexibility within a broad 
hydrophobic channel. Molecular dynamics simulations 
on the SM complex provide an explanation for the 
opposing impact these mutations had on hydrolysis of 
SM versus LPC. In particular, the distance of the sphin- 
goid tail to F80 increased during the simulation, 



stabilizing at a value over 8 A. This distance is suffi- 
ciently high that minimal interactions between SM and 
F80 occur, consistent with the catalytic activity similar 
to WT for the F80A mutant when tested using SM as 
the substrate. The SM sphingoid tail remained within 
contact distance of Y142 throughout the simulation 
(averaging 4.5 A), although the closest contact involved 
backbone atoms of Y142. The sidechain of Y142 forms 
contacts between strands of the p -sheet as well as to a 
neighboring helix. Mutation to a smaller alanine residue 
may therefore alter the shape or size of the binding 
pocket for the second tail of SM rather than removing a 
direct contact. The closest contact of the SM sphingoid 
tail averaged about 2.5 A with Y166 with a very low 
standard deviation of 0.2 A. A likely explanation for the 
enhanced activity of the Y166A mutant relative to WT 
for SM hydrolysis is that greater conformational flexibil- 
ity of the SM sphingoid tail would be accommodated by 
the change from a bulky tyrosine to a more compact 
alanine sidechain, reducing the entropic penalty of 
transferring SM into the NPP7 binding pocket. Figure 8 
shows a snapshot from the simulation emphasizing the 
close contact with Y166 and distance from F80. The 
rigid alkene portion of the sphingoid chain sits directly 
above the bulky sidechain of W119 (not shown to sim- 
plify view of mutated residue positions), preventing the 
sphingoid chain of SM from exercising the same flexibil- 
ity in the WT as the saturated sidechains of LPC 16:0 
and PAF 16:0. 

One mutation site, E169, was selected on the basis of 
isolated contact with the second tail of PAF (short 2:0 



Parrill et al. BMC Biochemistry 201 1, 12:65 
http://www.biomedcentral.com/1471-2091/12/65 



Page 12 of 15 




Figure 9 SM position in NPP7 relative to mutated residues. Structure shown represents an energy-minimized snapshot from the MD 
simulation. Distances from sidechains to SM are 2.4 A (Y166 to sphingoid tail), 2.5 A (L107 to choline), 2.9A (E169 to palmitoyl tail), 3.7 A (F275 
aromatic centroid to choline), 5.5 A (Y142 to sphingoid tail), and 6.9 A (F80 to sphingoid tail). Mutation of F141 (yellow) to serine was reported 
by Duan, et al. to substantially reduce sphingomyelinase activity of NPP7 [15]. Distance between F141 and the palmitoyl tail is 2.8 A. 
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chain) and SM (16:0 chain). E169 was mutated to both 
the similarly-sized glutamine (E169Q) and much smaller 
alanine (E169A) residue. The mutation to glutamine 
exhibited minimal differences from WT for LPC, PAF 
and pNPPC, but showed a 1.5-fold enhanced activity 
against SM. In contrast, the mutation to alanine showed 
substantially decreased activity for LPC hydrolysis and 
less drastic decreases for the remaining substrates. Fig- 
ure 8 shows that a hydrogen bond between the free gly- 
cerol hydroxyl group of LPC 16:0 and E169 forms 
during the molecular dynamics simulation of the LPC 
complex. This hydrogen bond could occur in the E169Q 
but not in the E169A mutant. Thus, the simulation 
results on the LPC complex with NPP7 are consistent 
with the observed decrease in activity only of the E169A 
mutant relative to WT. In contrast to the LPC complex, 
Figure 9 shows that the interaction of SM with E169 
predominantly involves contact between the fatty acyl 
chain and the E169 sidechain. Thus enhanced activity 
due to elimination of sidechain ionization in the E169Q 
mutant and relatively limited activity difference in the 
El 69 A mutant are consistent with the simulation results 
on the SM complex with NPP7. 

Discussion 

The choline headgroups of four different substrates are 
positioned consistently near F275 and L107 in prelimin- 
ary homology models of NPP7 (Figure 5). The L107F 
mutant was able to enhance the catalytic function of 
NPP7 against three of the four substrates tested, with 
reduced function only against the largest substrate due 
to steric crowding (Figure 7, and Tables 1, 2, 3, and 4). 
This finding was consistent with the preliminary models 
which suggested a cation-n interaction might form in 
the mutant to enhance interactions with the choline 
groups of the three smaller substrates (Figure 5). All 
other human NPP family members have a phenylalanine 
at the corresponding position. The observed reduction 
of SM hydrolysis upon mutation of L107 to phenylala- 
nine suggests that the presence of a phenylalanine in 
the other NPP family members contributes to substrate 
selectivity by steric interactions with the bulkier SM 
structure relative to other candidate substrates. This is 
unlikely to be the sole factor determining whether an 
NPP enzyme can utilize SM as a substrate as the L107F 
mutant did not completely lose the ability to hydrolyze 
SM. The second choline interaction site was probed 
using an F275A mutant, which showed impaired cataly- 
tic function against all four substrates as predicted by 
the preliminary models (Figure 5 and Tables 1, 2, 3, and 
4). This interaction site was previously described by 
Duan, et al, in their modeling study on NPP7. [15] They 
found that an F275G mutant showed almost no sphin- 
gomyelinase activity, consistent with our reduced 



activity of the F275A mutant. Perhaps surprisingly, an 
alignment of the mammalian NPP enzymes shows 
remarkable variability in the residue types that align 
with F275. Neither NPP2/ATX nor NPP6 have a pheny- 
lalanine aligned with F275, despite shared recognition of 
phospholipids with choline headgroups. NPP2/ATX 
exhibits D434 at the corresponding position, which may 
form a favorable ion pair with the cationic choline head- 
group instead of the cation-n interaction observed in 
NPP7. NPP6 aligns Q266 against F275 of NPP7. The 
role of Q266 in NPP6 will require further investigation, 
since a strongly favorable interaction with the choline 
headgroup is unlikely. Notably, both NPP4 and NPP5 
have tyrosine residues (260 and 261, respectively) corre- 
sponding to F275. Although the endogenous substrates 
of NPP4 and NPP5 are unknown, the tyrosine sidechain 
is capable of the same type of cation-n interaction as 
the phenylalanine sidechain. 

In contrast to the very high consistency between the 
headgroup position in the model and experimental 
mutagenesis results, the preliminary modeled positions 
of the substrate tails were not uniformly consistent with 
the experimental mutagenesis results. Molecular 
dynamics studies indicated that the hydrophobic tails 
showed dynamic interactions within a channel lined 
with hydrophobic residues. This dynamic behavior 
explained how mutation of F80 and Y166 to alanine 
could both have detrimental impact on LPC hydrolysis 
although they were unable to concurrently interact with 
the LPC hydrophobic tail (Figures 7 and 8). This finding 
compares well with the divergent positions of the LPA 
hydrophobic tails observed in the crystallographic struc- 
tures of NPP2/ATX (Figure 10) [6]. The hydrophobic 
tails of LPA species in PDB [18] entries 3NKQ[6] and 
3NKR [6] are located near Y214, which corresponds to 
F80 in NPP7. The hydrophobic tail of the LPA species 
in PDB entry 3NKP[6] contacts W275, which corre- 
sponds most closely in the superposition to NPP7 resi- 
due Y166 in the insertion loop first described by 
Nishimasu, et al. in their comparison of the bacterial 
NPP and NPP2 [6]. The hydrophobic tails of LPA spe- 
cies in PDB entries 3NKN[6] and 3NKO [6] are posi- 
tioned intermediately between these two sidechains. 

The experimental mutation results are compared in 
Figure 7, and clearly demonstrate a very different pat- 
tern for SM relative to LPC, PAF and pNPPC, which 
show more similar profiles. SM is the only substrate 
tested that has two long and flexible hydrophobic tails 
(Figure 1). In our model, these tails occupy completely 
different binding subpockets, above and below the con- 
served P -sheet core. This prediction is in direct contrast 
to the modeling studies of Duan, et al, in which both 
tails were placed into the region below the conserved P- 
sheet core [15]. In their homology modeling, the loop 
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Figure 10 Crystallographic positions of LPA species in NPP2/ATX (PDB entries 3NKN, 3NK0, 3NKP, 3NKQ, 3NKR[6]). For clarity, the 
backbone of only 3NKN is represented as a ribbon colored by secondary structure (yellow for (3-sheet, red for a-helix, and blue for turns) and 
the zinc ions of only 3NKN are represented as cyan-colored spheres. Residues Y214 and W275 are shown as sticks due to the correspondence 
with F80 and Y166 of NPP7, respectively. 



segment from F140 to E171 was subjected to refinement 
to generate an ensemble of loop conformations. Confor- 
mations were selected for further modeling that kept 
this lower binding subpocket more open to allow both 
chains to fit together. Our mutation of E169Q, which 
showed 1.5-fold enhanced hydrolysis of SM with rela- 
tively little impact on hydrolysis of the remaining sub- 
strates, suggests that the second tail of SM is not 
occupying the same site as the hydrophobic tails of the 
single-tail substrates, as more similarity of mutagenesis 
results between SM and other substrates would be 
expected. In fact, we interpret the enhanced activity of 
the E169Q mutant against SM as a reduction in the 
charge allowing improved hydrophobic interactions with 
the second SM tail. Furthermore, Figure 9 shows that 
the palmitoyl tail of SM makes close van derWaals con- 
tact with the sidechain of F141, consistent with poor 
sphingomyelinase activity observed for an F141S mutant 
by Duan, et al, a finding that was not anticipated based 
on the common positioning of the hydrophobic tails of 
SM in their model below the conserved p -sheet core 
[15]. Based on the current results, the F141S mutant is 
likely to have activity similar to WT NPP7 against sin- 
gle-tail substrates such as LPC. Importantly, as SM is 
the optimal substrate for NPP7 out of those tested to 
date, it is inconsistent to expect that the conformational 
flexibility of the nonpolar tails would be restricted by 



common occupation of a pocket of limited size due to 
the unfavorable entropic contribution that the loss of 
flexibility would produce. Thus several residues in the 
longer insertion loop found in NPP7 relative to other 
mammalian and bacterial NPP family members have 
been demonstrated to play a role in substrate recogni- 
tion, and likely contribute to substrate selectivity as pro- 
posed by Nishimasu, et al. [6] based on their 
comparison of bacterial NPP and mouse NPP2. 

Conclusions 

The current study has defined important differences in 
recognition of SM and other phospholipid substrates of 
NPP7. Two contributions to the unique activity of NPP7 
against SM have been identified. First is the importance 
of the amino acid at position 107. NPP7 has leucine at 
this position although related NPP enzymes have pheny- 
lalanine at the corresponding positions. The detrimental 
impact of the NPP7-L107F mutant on hydrolysis of only 
SM, compared to the positive impacts on the hydrolysis 
of LPC, PAF and pNPPC supports our conclusion that 
this position contributes to substrate selectivity within 
the NPP family. The substrate-specific impacts of the 
mutations at position E169 indicate a second contribu- 
tion to the unique recognition of SM, a second hydro- 
phobic binding pocket for the second hydrophobic tail 
that is unique to SM. 
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Additional material 



Additional file 1: Substrate complexes shown in Figure 3. Atomic 
coordinates of selected homology model and initial docked complexes 
of four substrates. 



List of abbreviations 

ATX: Autotaxin; CCM: Concentrated conditioned media; DMEM: Dulbecco's 
modified eagle's medium; HEK: Human embryonic kidney; LPC: 
Lysophosphatidylcholine; LPL: Lysophospholipase; MD: Molecular dynamics; 
MOE: Molecular Operating Environment; MWCO: Molecular weight cutoff; 
NFDM: Non-fat dry milk; NPP: Nucleotide pyrophosphatase/ 
phosphodiesterase; PAF: Platelet activating factor; PCR: Polymerase chain 
reaction; PDB: Protein Data Bank; pNPPC: para-nitrophenylphosphorylcholine; 
PVDM: polyvinylidine difluoride; RMSD: Root mean square deviation; RMSG: 
Root mean square gradient; SDS-PAGE: Sodium dodecyl sulfate 
polyacrylamide gel electrophoresis; SM: sphingomyelin; SPC: 
sphingosylphosphorylcholine; TBST: Tris-buffered saline with TWEEN. 

Acknowledgements and Funding 

A mammalian expression vector pcDNA4/T0/myc His B containing the full- 
length human NPP7 sequence was a generous gift from Dr. Rui-Dong Duan 
(Lund University, Sweden). The MOE software was generously provided by 
the Chemical Computing Group (Montreal, Canada). 

Author details 

department of Chemistry and the Computational Research on Materials 
Institute, The University of Memphis, Memphis, TN 38152, USA. department 
of Chemistry, The University of Memphis, Memphis, TN 38152, USA. 

Authors' contributions 

ALP performed the molecular dynamics studies, participated in the 
development of homology models, contributed to the conception of the 
study, critically reviewed the analysis of kinetics assays, and completed the 
manuscript. IWJ participated in the conception of the study, participated in 
the development of homology models and substrate complexes, performed 
the kinetic assays, analyzed the kinetic assay results, and wrote an initial 
partial draft of the manuscript. TCTP participated in the conception of the 
study, performed the control studies shown in Figure 6, and critically 
reviewed the analysis of kinetics assay results. DLB participated in the 
conception of the study and critically reviewed all results. All authors read 
and approved the final manuscript. 

Received: 14 July 201 1 Accepted: 16 December 2011 
Published: 16 December 2011 

References 

1. Gijsbers R, Aoki J, Arai H, Bollen M: The hydrolysis of lysophospholipids 
and nucleotides by autotaxin (NPP2) involves a single catalytic site. FEBS 
Lett 2003, 538(1 -3):60-64. 

2. Duan RD, Bergman T, Xu N, Wu J, Cheng Y, Duan J, Nelander S, 
Palmberg C, Nilsson A: Identification of human intestinal alkaline 
sphingomyelinase as a novel ecto-enzyme related to the nucleotide 
phosphodiesterase family. J Biol Chem 2003, 278(40):38528-38536. 

3. Sakagami H, Aoki J, Natori Y, Nishikawa K, Kakehi Y, Arai H: Biochemical 
and molecular characterization of a novel choline-specific 
glycerophosphodiester phosphodiesterase belonging to the nucleotide 
pyrophosphatase/phosphodiesterase (NPP) family. J Biol Chem 2005. 

4. Stefan C, Jansen S, Bollen M: NPP-type ectophosphodiesterases: unity in 
diversity. Trends Biochem Sci 2005, 30(10)542-550. 

5. Hausmann J, Kamtekar S, Christodoulou E, Day JE, Wu T, Fulkerson Z, 
Albers HM, van Meeteren LA, Houben AJ, van Zeijl L, et al: Structural basis 
of substrate discrimination and integrin binding by autotaxin. Nat Struct 
Mol Biol 2011, 18(2):1 98-204. 

6. Nishimasu H, Okudaira S, Hama K, Mihara E, Dohmae N, Inoue A, Ishitani R, 
Takagi J, Aoki J, Nureki O: Crystal structure of autotaxin and insight into 
GPCR activation by lipid mediators. Nat Struct Mol Biol 201 1, 
18(2)205-212. 



7. Bollen M, Gijsbers R, Ceulemans H, Stalmans W, Stefan C: Nucleotide 
pyrophosphatases/phosphodiesterases on the move. Crit Rev Biochem 
Mol Biol 2000, 35(6):393-432. 

8. Goding JW, Grobben B, Siegers H: Physiological and pathophysiological 
functions of the ecto-nucleotide pyrophosphatase/phosphodiesterase 
family. Biochim Biophys Acta 2003, 1 638(1 ):1 -19. 

9. Cheng Y, Wu J, Hertervig E, Lindgren S, Duan D, Nilsson A, Duan RD: 
Identification of aberrant forms of alkaline sphingomyelinase (NPP7) 
associated with human liver tumorigenesis. Br J Cancer 2007, 

97(1 0):1 441 -1448. 

10. Wu J, Nilsson A, Jonsson BA, Stenstad H, Agace W, Cheng Y, Duan RD: 
Intestinal alkaline sphingomyelinase hydrolyses and inactivates platelet- 
activating factor by a phospholipase C activity. Biochem J 2006, 394(Pt 
1):299-308. 

11. Duan RD, Nilsson A: Metabolism of sphingolipids in the gut and its 
relation to inflammation and cancer development. Prog Lipid Res 2009, 

48(1):62-72. 

12. Zhang Y, Cheng Y, Hansen GH, Niels-Christiansen LL, Koentgen F, 
Ohlsson L, Nilsson A, Duan RD: Crucial role of alkaline sphingomyelinase 
in sphingomyelin digestion: a study on enzyme knockout mice. J Lipid 
Res 2011, 52(4)771-781. 

13. Wu J, Cheng Y, Nilsson A, Duan RD: Identification of one exon deletion of 
intestinal alkaline sphingomyelinase in colon cancer HT-29 cells and a 
differentiation-related expression of the wild-type enzyme in Caco-2 
cells. Carcinogenesis 2004, 25(8):1 327-1 333. 

14. Wu J, Hansen GH, Nilsson A, Duan RD: Functional studies of human 
intestinal alkaline sphingomyelinase by deglycosylation and 
mutagenesis. Biochem J 2005, 386(Pt 1 ):1 53-1 60. 

15. Duan J, Wu J, Cheng Y, Duan RD: Understanding the molecular activity of 
alkaline sphingomyelinase (NPP7) by computer modeling. Biochemistry 
2010, 49(42):9096-9105. 

16. Zalatan JG, Fenn TD, Brunger AT, Herschlag D: Structural and functional 
comparisons of nucleotide pyrophosphatase/phosphodiesterase and 
alkaline phosphatase: implications for mechanism and evolution. 
Biochemistry 2006, 45(32):9788-9803. 

17. MOE: Montreal: Chemical Computing Group;, 2009.10 2010. 

18. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, 
Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res 2000, 
28:235-242. 

19. Halgren TA: Merck Molecular Force Field. I. Basis, Form, Scope, 
Parameterization, and Performance of MMFF94*. J Comput Chem 1996, 
17(5 & 6):490-519. 

20. Onufriev A, Bashford D, Case DA: Modification of the Generalized Born 
Model Suitable for Macromolecules. J Phys Chem B 2000, 

104(15)3712-3720. 

21. Bond SD, Leimkuhler BJ, Laird BB: The Nose-Poincare Method for Constant 
Temperature Molecular Dynamics. J Comput Phys 1999, 151(1 ):1 14-134. 

22. Sturgeon JB, Laird BB: Symplectic algorithm for constant-pressure 
molecular dynamics using a Nose-Poincare thermostate. J Chem Phys 
2000, 1 1 2(8)3474-3482. 

23. Wu J, Liu F, Nilsson A, Duan RD: Pancreatic trypsin cleaves intestinal 
alkaline sphingomyelinase from mucosa and enhances the 
sphingomyelinase activity. Am J Physiol Gastrointest Liver Physiol 2004, 
287(5):G967-973. 

24. Suresh PS, Olubiyi O, Thirunavukkarasu C, Strodel B, Kumar MS: Molecular 
modeling of human alkaline sphingomyelinase. Bioinformation 2011, 
6(2):78-82. 



doi:1 0.1 1 86/1 471 -2091 -1 2-65 

Cite this article as: Parrill et al:. Computational identification and 
experimental characterization of substrate binding determinants of 
nucleotide pyrophosphatase/phosphodiesterase 7. BMC Biochemistry 

2011 12:65. 



